۰ 0۳0۵ ۲۱۵۱0 ۵ امصسهمل صهتصره] 
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+52۱6 ۷۲۰ ۳ )صحطهل ۱۲۱ و" تاماافطه]-تعه۱۲ ۷۲۰ ۳ نافططعظ ٩:‏ 


:۰ و۱ ان 0 ۳0۷ | . 27-01-2022 :۲26661۷60 

لمءزومامممط همصننام ,(۱۷۲,,)2023 رتطملهی ی ,.۱۷ مصفطعل ,۷ بتاهالمطه]۷-توعهاط بیط متلهاطاع۴ 21-05-02 :۲6۷1860 
,(1) 21 ,6560۲6 ۳۵0۵6 10ع۳۱ ۵۴ ۵۵۳۵( ۱۲۵۳/۵۵۵۲ بعاجنا ۳۱۵۱۵ و۵ 2860ظ عصتام۱۷]۵0 )صعصوم م1۵۷۵ 19-06-22 :۸6660۱60 
5 ۵ .()0۳20 205 طفنام‌صظ ط)ز< صوتوعظ ص) .47-60 


ار ۱۱۱۹ 


م2 2تصناه؟ مصلودع6۲ض1 ر‌عصمهافصا ۳۵۶ ممتامن0متم وم ۵۶ فصمتانلومم فطا عصتلهعع0 ب«النمه: قوذ ععصقطه ماقصتت 
4 وم0ع۵و )ههام -عمفعتاو 89 رعمممن008860 2 ۸۵ ,۷۷۵۲۱0 معط ۵۶ قاتوج ]وم ما مقهعععض1 ما 101۲6628000 و1 اند 230 
۷ 2۳ ]00 اوه ملعم اصهتعاه)-ووعتاه .عتلنام‌تنيه معتابظ ۶۵۲ ۵0وی قصح ‏ 106201160 ۵ )فسظ وعمرامصمع 
معط ۵۶ توحاصصمجهر ۵ و1 ,۷۷۱۱۱۵ همصنتان حصتلمممصمفط وه موم ماه ردمصتنام .0)60ع1قفص مضه 11260 ارع0صا 
عص ما الم ۲۵00 مهن له همصتیان گم حمتاتتابه 4صح صمته فان فطا مصتامطصمظ .اتحصه مهععمط)صه‌تمصض 
0 .86608 1000 و0تصصدرمی فظ) ۵۶ )ت22 ۵۲0۷106 0ج ومحطرمعصا فتعصته؟ مفهعتمصا رصمتامت0مم فاطاقصته) فتاه روتصامم 
ع عصتصتماع رنه مفهتم200۲0 معط عمتاهعامو فج میاه ومتلباه فیامتیه۷ 107 ۱۷960 مععه 2۷۵ قام0معظ ممتاهلباظله 
۲ ۲60 طممعوع۱ ۲۱۵۱0 ,ظ)«مع وم معصقطل متقصصتاه مه پاتقعنل ۵۶ ]0ع1ه مطا مصتاهت0عتم رماع عصتاصهان اوعط 
۵ هه اه تتمحرمحر. صح فصصت 92۷6 صقه فاعم ممتاهانته. عمجم فلتطاه رتجمحمصهر. مضه مصصتاً 0۶ 
۰ صنه ل2اصمصتتممرم 


۱۷]۲۵۸)۵۳1215 2:0 ۵۹ 


2 ۶ 10 عطا ما عون مامتهومز 10 ۳1۵ 0۲۵۷1۵۵۵ ۷۵2۵0 ۵۶ وصمنعع۲ ۵ صا 0ع۲عنصمن ۱728 طمت2عوع2 فقط 1" 
وعصنا 011160مهصه مصتوتصمتم 5 1010060 وعمامه؟ له)صمصته‌مب .فممتاهه‌نامم۲ 3 ط۳ صعتعع0 »0۱۵01 میعاممصمع ۲۵00۴2۲26۵0 
,6ص زاتتامحط م0126مصصعاصا عنام ۵ 60افلدمم ومصاً م1 .هنن همهع۲10 از تعتصمن طمتهمعم بمتصتلمه ۷220 ظ1 
4 ۲۵0وصاصصیه مصنا تمد متقا عم ررم۸من 0۵۲9۴ 6 220 (رتامن ۲9 3 ]0۱9 2 ,رتم۲9۳۴ 1 ۲۵0هواصناط 
0 ۷۵۲ عصلله]-عامظ مضه عصتام‌صهه .5 0عتمطصصه هم۲1062. مد تاه اتستمطظ «انده فط فصه ,رظن 9۴ظ) 
1 ۳ .1 جموه ۵۶ قمعهاه لهءنعمامممطن مطا ۵۶ وومتو0۳۵ مطا ما ممت0۲00۵0۲ طا روججل عععصا وه ومصم رزاتقلناعع۲ 
۶ 02۷7 00۲۲ عط) م) مبال ,عاهد) طاعصع1 عمط هصتنهمهعم م1 .عمهناع‌ص12 ۲1 ها ۵۲۵۵۲0 ق۵ ظ)۲66-02۲-2۲0۷9عع0 هم 2860 
6 موتتامط 13.6 تفه 200 رقصم م۷2 م۲2 اممصه۷۵۱0۵ع0 فطا روعتامط 12.5 هه هبا معط فط طد ففصتاً له 10۲ رهممتنم 
فصح اهتنا وه« [06مظ فطع مرصفط1. 0 2 و772 [ع0مصظ عط صد تمصع موجه م1 .2۵۲0 قه۷ ماه کصمصصم‌ماع0۵۷ 
۰ 3۱ مگ صم‌لم) 0262 طاا مهبم 


«مع) 6( 200 و1666 


4 10 4.4 60۷۵۵۵ (۲ا۴۳]9) عم متفتاوو قمع ۲۵۵۲ و126 1۵ 7 ۵07۵۵ ۷۵/۵0 ۵۶ اممتمتاههه (0۱۷) ۲۷۲9۳۲ 

محعصصر مطا ,0.98 10 0.96 ومع9۵07 (ظ) پمممتم‌تلله ۵01 ر 4۵ ۵0.99 ۵۵07۵۵8 (0) 146 اصمصصمم‌تعه ۷۷1180 ,025 
0 ۲۷۷۵۵ 09 ماع 0۲ امملماااه0 مط] 20 0.08 6۵ 0.05 ممع6007ه قه/۱۷ (ظ۲۷۱) ۵61ص فط مصمن ممتا0۵۷12 
ص حمتت 061محظ مطا 16 همصنتان ۵۶ 11۱07 ما چم قطا ۵۶ متقصصتایی 2000 2 1۳0162060 ومباله۷ 1۳8686" .9846 0 66 92 
۶ ]0۵1110160 1۳ .۷۵11 تقووعععظ مط) 0مصنهع 0ماهناصله مصله1102 ما سمل ۶ وعبالم۷ فط هه رععهتاوصها ۳9 
7 2004 6.2 0۵0۷۵۵۵ و1 (رفا ۵ مه متمناوو صقفحظط )۲۵0 مطا رم٩8.6‏ م۵ 6.8 ممما۵ه ود (0۲) طرفاکه ۵ صمتامه۱۷ 
۷۵۵6 (ظ۷۲) ۵061 6 مطم جمتاه06۷1 عمط م1 ,0.92 220 0.75 0۵0۳7۵۵ و1 (0) «ع10 موه )مصصلز ۷۷ معط ر0215 
205 ۲0۵86 .9896 280 92۵ 06۷۵۵8 ۷۷۵6 09 ممتحعصتصمعل ۵ اصم‌تلههه فطا 20 0.08 1۵0 0.05 صمع۵07 
فص ,موماع1۱۸۵ ۲91 عطا ص صع)۳/۲1 06مص فط طز تاتتصفظ لممتعماملو ردان م) چدل فطا ۵۶ متقصتاوه 2000 2 1۳01621060 


,۱۷۱2۵20 ,مططوه ۶ تهب نوت تدهم0ع همه ۵۶ امد ,برع مامصطهعاهته۸ ۵۶ اممصاندمعنا رمتقتمته .بط -1 
1122 

مد ,لمططوه]۱۷ رمططفه ۱ ۵ نمهب نما لد 00ع۱ تیوه که تمد ررعمامصطهعاهتوظ ۵۶ ]ممحصا 12۳6۵2۲ ,۳۲۵0169501 -2 
۵ ۲۵85108 مه ممتمع 0 بطم‌تههومک لمعبهانم‌تبوه رتماجعن امعم اتصتلمه لمممتاداط ,۲۳۵16950۲ )صهاوزو۸۵ -3 
,۷۲۵70 ,(۲۵ض) 
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۸ نشریه پژوهشهای زراعی ایران. جلد ۲۱ شماره ۱ بهار ۱2۰۲ 


ممتامله ۵۱ مه صمتاهتطانلمی .ت۷۵ «تدووعععه مطا 4مصنمع باتتاممط آمعتومامتووطم 0مماناصته فطا مه فمتتله۷ رهق فطا 
(0۷) مماهه۵ن و۴ ۵ ممتامنه۷ ۵۴ هه ,رتاوا۵؟) مه معمننو صوعجط ۲۵0۵۲ مصتفدا پ«مصماهاله عم 0 
0 0 )0010160 ,(ظ۷۲) 06۷1200 مه صهممه ره ورطظا) «مصمته‌تاه م۵0 ورق) 106 اصممطرعمع2 ۱۷۷۱۱82۵۲ 

۰ اوه مه تافص آمءنعمآمتورطم 0ممع فمه عصنه/0 بصمتامصتصهع آتاصه رم 1 :1 ]معا مصنا بک) 


صمتعصل‌صهم) 


062۲66۵-02۷۲ ۵۶ فصصعا ما مصننان 10۲ 0۲0۵۵2۲60 ۵061 مصلنان مطا بتقطا 10162000 ناو فنط ۵۶ فلناومر فط 1 
تسه 0۶ عصنعها ما کصقام فنطا ۵۶ رحمتاح‌ستهمه هه مصته۳ 110 ,ععممع‌مصصهم) وعهماه لماممصمماه7ع0 مطا ماهنلع01 ۷۷۵11 
ون .عوعته )صمعللنل صا مق مصتاصهام مامتيممتممه عط 0عصتصیعاع0 مامط و 0 صق مضه (عفق1 هه صتتعص ررآتقع) 
10۴ 011600 ۵100ممامصه 0ص ممنهمموصها امومع ماقتاله7ه 6۵ ۵0وبا 9 07 صقه 061مهصه-ماناه 0ماهتصانلهه 
متقصناه معط 0ص اصمسه صا فصماوروه ممتاد تایه همصتیین طا فصو طامع ۵۶ موصه 10 2 طا عصل۳22۱ 
0۰ ۱) ما ۱۷۵ م0عناممه مه طمدعوم- م6020 ها 60وی وه و ۵061ج-ماناو فتطا رم۵۲6۴0ظ1 .فصمتاتصمم 


مق م0۳ ,ودوج معتوع0 م00۳1 ,۲9 :1۵۲۱۷۵۲۵۵ 


نشربه پژوهشهای زراعی ایران 
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حلد ۰۱ شماره 5 بهار ۲ص ۶۷-۰ 


شبیه‌سازی نمو فنولوژیک کینوا بر اساس داده‌های مزرعه‌ای 


شهاب اقبالی شاه آیاد۵ ال مهدی نصیری محلاتی 69 محسن جهان ۵ معصو مه صالحی ۲ 


تاریخ دریافت: ۱۴۰۰/۱۱/۰۷ 


تاریخ پذیرش: ۱۴۰۱/۰۳/۲۹ 


چکیده 


هدف از این پژوهش تهیه مدلی برای شبیه‌سازی نمو کینوا (.۷۷1110 7100 6۱0۴(00:10:/)) و واسنجی و تعیین اعتبار آن بر اساس داده‌های 
مزرعه‌ای است. این تحقیق در دو منطقه استان یزد با ۱۰ آزمایش جداگانه و در قالب طرح بلوک‌های کامل تصادفی با ۳ تکرار انجام شد. عامل‌های 


آزمایش شامل ۵ لاین امید بخش اصلاح‌شده در مرکز تحقیقات شوری یزد به‌همراه رقم تی‌تی‌کاکا بود. نمونه‌گیری و یادداشت‌برداری‌ها به‌طور مرتب» 


هر سه روز یک‌بار و متناسب با پیشرفت مراحل فنولوژیک هر لاین انجام شد. مدلی بر اساس درجه-روز-رشد با زبان برنامه‌نویسی 91 تهیه و سپس 
اقدام به واسنجی و ارزیابی مدل با داده‌های برداشت شده از مزرعه شد. واسنجی و ارزیابی کارایی مدل با استفاده از شاخص‌های جذر میانگین مربعات 
خطا (1116۳5» ضریب تغییرات یا درصد 01271517 ( 46۳۷ شاخص توافق ویلموت (0)» کارایی مدل (0۷13 میانگین انحراف از مدل (418 ضریب 
تبیین (**4 آزمون خط ۱:۱ روز تا سبزشدن, گل‌دهی و رسیدگی فیزیولوژیک خوب برآورد شد. این مدل تعبین اعتبارشده و هم‌اکنون می‌تواند برای 
ارزیابی تاثیرات مختلف دمایی و فتوپریودی برای تصمیم‌گیری در طیف وسیعی از محیط‌های رشد در نظام‌های کشت کینوا در شرایط اقلیمی جاری و 
اینده مورد استفاده قرار گیرد. بنابراین از این مدل در کارهای آموزشی- تحقیقاتی و کاربردی در مزرعه می‌توان استفاده کرد. 


واژه‌های کلیدی: تاریخ کاشت. ۳9 درجه-روز-رشد 


مقدمه 

کینوا (۷۷1110 04 ۵70۴۵0۵:)) گیاهی از خانواده 
۵ است که بومی کوه‌های آند در بولیوی» شیلی و 
پرو بوده و در حال حاضر بولیوی» پرو و اکوادور مهم‌ترین 
تولیدکنندگان اين گیاه هستند (2010 ,. 6 021۷762 - ۷۵۵2). این 
گیاه انعطاف‌پذیری بالایی در سازگاری به فتوپریود. ارتفاع از سطح 
دریاء دما و خصوصیات خاک دارد؛ به‌طوری که» در مناطقی با ارتفاع 
۰-۰ متر از سطح دریاء 1 بین ۶ تا ۸/۵ و مناطقی با اقلیم 
خشک» نیمه‌خشک» حاره‌ای و مرطوب قابلیت کشت دارد. دوره رشد 
کینوا بسته به منطقه و نوع رقم مورد کشت بین ۷۰ تا ۲۴۰ روز متفیر 


۱- دانشآموخته دکتری اگرواکولوژی, دانشکده کشاورزی دانشگاه فردوسی مشهد. 

مشهد ایران 

۲- استاد گروه اگرواکولوژی دانشکده کشاورزی» دانشگاه فردوسی مشهد. مشهد 

ابا 

۳- استادیار مرکز ملی تحقیقات شوری» سازمان تحقیقات آموزش و ترویج 

کشاورزی» یزده ایران 

(*:- نویسنده مسئول: :۱۷۵۵00۵۱۷۵6۵۵0 _داه :اتفصصط) 
5 /005://101.012/10.22067 ۱ 


است (2016 ,.۵1 6 ۷22). محصول اصلی این گیاه دانه است که 
ارزش غذایی بالایی ازنظر پروتئین» توازن اسیدآمینه» چربی غیراشباع, 
ویتامین‌ها و مواد معدنی دارد (2017 ,4 2 ۳1506۲). بنابر 
گزارش‌های سازمان خواروبار جهانی این گیاه حتی در اراضی فقیر و 
غیرقابل کشت نیز می‌توان کشت کرد. علاوه بر لین کشت کینوا در 
تلوب زراعی با غلات به منظور کنترل علف‌های هرز و بیماری‌ها 
توصیه می‌شود. ترویج کشت و تغذیه کینوا باعث تنوع‌بخشی 
محصولات غذایی در کشور تولید پایدار افزایش درآمد کشاورزان و 
تأمین بخشی از نیاز غذای جامعه خواهد شد. با توجه به نقش بالقوه 
این گیاه در بهبود امنیت غذایی در آینده و تحمل بالا به تنش‌های 
شوری و خشکی سازمان خوار و بار جهانی سال ۲۰۱۳ را به نام این 
گیاه نام گذاری کرد (2014 ,۳۸0 

مدل ریاضیی معادله یا مجموعه‌ای از معادلات است که رفتار هر 
سیستم را به‌طور کمی توصیف می‌کند. سیستم, به بخش محدودی از 
عقاه تاقنی اطاات مش ان ‌های اس ی ات آه 
متقابل با هم داشته و در طی زمان تغییر می‌کند. از مدل‌های 
شبیه‌سازی گیاهان زراعی برای انحام مطالعات مختلف از جمله انتخاب 
گیاه و رقم مناسب تعیین خصوصیات مطلوب گیاهی تعیین بهترین 


۵۰ نشریه پژوهشهای زراعی ایران. جلد ۱ شماره (۵ ببهار رف 


مدیریت زراعی» پیش‌بینی اثر تنوع و تغییر اقلیم بر رشد» تخمین 
عملکرد. برآورد ظرفیت تولید منطقه‌ای» تعیین خطمشی برای 
به‌نژادی» تعیین اولویت‌های تحقیقاتی و انتقال تکنولوژی طبقه‌بندی 
آگرواکولوژیک استفاده شده است. اجرای تحقیقات مزرعه‌ای به صرف 
وقت و هزینه زیاد نیاز دارد در حالی‌که مدل‌های شبیه‌سازی رایانه‌ای 
می‌تواند با انجام شبیه‌سازی آزمایش‌های وسیع باعث صرفه‌جویی در 
وقت و هزینه شوند (2013 ,.۵1 21 ۷۲2006۷۷5 اما باید توجه داشت 
که مدل‌های شبیه‌سازی گیاهان زراعی به هیچ‌وجه نمی‌توانند 
جایگزین آزمایشات مزرعه‌ای شوند بلکه بهترین حالت آن است که 
به‌عنوان مکمل این آزمایش‌ها استفاده شوند. در این صورت مدل‌های 
گیاهان زراعی خواهند توانست تا حدودی نیاز به آزمایشات مزرعه را 
کاهش دهند و اثربخشی نتایج این آزمایش‌هارا بیشتر کنند. 
مدل‌سازی در آموزش‌های کاربردی برای محققین, دانشجویان» 
کارشناسان کشاورزی و حتی بهره‌برداران. ضمن کاهش هزینه‌ها 
موجب شناخت و درک بهتر استفاده‌کنندگان از فرایندهای پبجیده 
زیستی, فیزیولوژی و اکولوژی گیاهان زراعی خواهد شد؛ تابه حال 
مدل‌های متفاوتی برای پیش‌بینی فنولوژی گیاهان زراعی مختلف 
استفاده شده است ,4 ۶ ۲6۲/6۲0 :2021 ,اه 6 طدحاطل) 
(2006 ,.اه ۶۶ تصم‌آه5 :1999 

تاریخ کاشت» نسبت به سایر روش‌های مدیریت زراعی بیشترین 
تأثیر را بر ویژگی‌های فنولوژیکی و رشدی گیاهان می‌گذارد 
(2018 ,اه ۵1 صدنصعدهطع۳0۳) و نقش مهمی در بهبود درصد 
سبز و رسیدن به حداکثر عملکرد دارد. تاریخ کشت مناسب کینوا به 
اقلیم منطقه و رقم بستگی دارد. ارقام زودرس کینوا برای کشت در 
شهریور و پا کشت بهاره. ارقام متوسطرس کینوا برای کشت در مرداد 
با کیت کارا کم فیررسن کشا مرای کس در اسف با 
کف اسان در ستاطق مد ششک کسور مایب هت که 
(2021 ,./۸ 6 تل02ع۳). این ارقام در مناطق معتدله در صورتی که 
در مرحله ۱۳ برگی باشند در طول زمستان زنده می‌مانند ولی بذر 
تولید نمی کنند (2017 ,نصهداعع(1 ک ,طع621). تعیین مناسب‌ترین 
تاریخ کاشت یعنی تعیین زمانی که رشد رویشی و زایشی گیاه بیشترین 
هماهنگی را با عوامل اقلیمی داشته و گیاه کمتر با شرایط نامساعد 
محیطی روبه‌رو شود (2003 ,.۵1 61 0۵۱011007). 

پیش‌بینی صحیح مراحل فنولوژی گیاهان زراعی برای بهینه‌سازی 
فعالیت‌های مدیریتی در مزرعه و سازگاری بهتر تقویم زراعی با 
سیست‌های آگرواکولتژیک خاضص بسیار میم سی‌باشد ( زم #ع 
4 0 جوان‌زنی و سبزشدن از مهم‌ترین فازهای فنولوژیکی 
استقرار گیاه است. جوانه‌زنی بذر فرایند پیجیده‌ای از فرایندها و 
واکنش‌ها می‌باشد که همگی تحت تاثیر دما انجام می‌شود. مهم‌ترین 
مرحله فنولوژیک, مرحله گل‌دهی یعنی عبور گیاه از مرحله رویشی به 
زایشی است. با توجه به این که بسیاری از فرایندهای فیزیولوژیک و 


مورفولوژیک با نوع گیاه تغییر می‌کند» کمی کردن نمو فنولوژیک در 
کلیه مدل‌های شبیه‌سازی گیاهان زراعی لازم است ( .01 6 )501 
06 گیاه کینوا بعد از گرده‌افشانی» وارد فاز شیری یا خمیری و 
رسیدگی‌دانه می‌شود. با پرشدن دانه کاسبرگ‌ها باز می‌شوند و 
به‌صورت ستاره‌ای درمی‌آیند و رنگ پوسته بیرونی بذر مشخص 
مي شود 

برخی محققین مدل‌سازی پاسخ گل‌دهی ٩‏ رقم کینوا به دماو 
فتوپریود را با رگرسیون خطی انجام دادند. ابتدا اثر دما و فتوپربود را 
جداگانه مورد بررسی قرار دادند و فتوپریود آستانه. فتوپریود بحرانی» 
دمای پایه و دمای بهینه گل‌دهی را شبیه‌سازی کردند و سپس از 
پارامترهای به‌دست آمده برای مدل‌سازی اثر متقابل دما و فتوپریود 
استفاده کردند. آن‌ها گزارش کردند بدون اثر متقابل دما و فتوپریود 
زمانی که دما بالاتر از ۲۰ درجه سانتی‌گراد بود مدل نمی‌توانست زمان 
گل‌دهی را خوب پیش‌بینی نماید. همچنین بدون اثر متقابل در دوره 
رویشی ارقام کینو؛ فتوپریود با عرض جغرافیایی منشا رقم رابطه منفی 
و با مدت زمان رشد رویشی رقم رابطه مثبت داشت. در بررسی اثر 
متقابل پارامترهای پاسخ دما و فتوپریود با عرض جغرافیایی منشا رقم 
اثر معنی‌دار نداشت (19990 ,۲1070عظ). 

نتایج آزمایشی در کرج نشان داد که کلیه ژنوتیپ‌های کینوا مورد 
ارزیابی برای کشت بهاره در این منطقه سازگاری داشتند. کشت در 
نیمه دوم اردیبهشت سال ۱۳۹۷ باعث مصادف شدن مرحله 
گرده‌افشانی و پرشدن دانه با دمای بالای هوا در تیرماه و مردادماه شده 
و باعث کاهش عملکرد دانه و اجزای عملکرد شد. رسیدگی‌دانه‌های 
کینوا در ماه‌های مرداد و شهریور, به دلیل کمبود غذای کافی برای 
پرندگان در این مقطع زمانی» خسارت هجوم پرندگان را در پی داشست 
(2018 ,۵ ۶ تحدطش). 

مدل فنولوژی» زمان وقوع یک فرآیند خاص در طی نمو 
پیش‌بینی می‌کند. به منظور مدیریت آفات تنظیم برنامه آبیاری و 
برخی دیگر از مسائل دوران کاشت گیاه لازم است سرعت رشد و نمو 
گیاه تعیین شود. سرعت نمو گیاهان زراعی از زمان کاشت تا زمان 
رسیدگی» شدیداً وابسته به دمای محیطی است که این گیاهان با آن 
مواجه می‌شوند. ازآنجا که شرایط آب و هوایی در طول سال‌های 
مختلف بسیار متنوع است بنابراین تاریخ‌های تقویمی نمی‌توانند 
مبنای خوبی جهت تصمیم گیری‌های مدیریتی در رابطه با موجودات 
و مراحل رشد و نمو آن‌ها باشند. بدین ترتیب با اندازه‌گیری مقدار 
گرمای تجمع‌یافته در طی زمان» یک قیاس زمانی فیزیولوژیک 
به‌دست می‌آید که ازنظر معیارهای زیست‌شناختی از دقت بیشتری 
نسبت به تاریخ تقویمی برخوردار است و این معیار در گياهان مختلف 
به نوع گیاه بستگی داشته و در ژنوتیپ‌های مختلف متفاوت است. 
برای محاسبه واحدهای گرمایی از روش درجه-روز-رشد یا ما6 
استفاده می‌شود. درجه -روز-رشد مجموع میانگین درجه حرارت بالاتر 
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از دمای پایه است. 

زمان حرارتی مورد نیاز برای تکمیل نمو کینوا در تاریخ‌های 
کاشت مختلف بین ۲۶۶۲ تا ۱۸۹۳ درجه روز گزارش کردند. دامنه 
تغییرات این زمان برای دوره سبز شدن بین ۸۲ تا ۱۵۰ مرحله 
غنچه‌دهی بین ۴۷۷ تا ۶۵۶ مرحله گرده‌افشانی ۳۴۸ تا ۸٩۰‏ و مرحله 
رسیدگی فیزیولوژیک ۱۴۳ تا ۱۱۲۱ درجه روز برآورد شد ( 61 نطهله6 
9 ,.1م). 

مدل‌های شبیه‌سازی نمو در دو گروه هستند» برخی مراحل نمو را 
می‌کنند و گروه دیگر پیشرفت نمو را بر اساس درجه روزهای تجمعی 
پیش‌بینی می‌کنند (2000 .۸7 ۵ 255171)؛ که مدل تهیه شده در 
این تحقیق از نوع دوم است. شبیه‌سازی نمو بخش مهمی از مدل‌های 
شبیه‌سازی رشد گیاهان از جمله کینوا است و تهیه» واسنجی و 
اعتبارسنجی چنین زیر مدل‌هایی برای ساخت مدل کامل شبیه‌سازی 
رشد کینوا ضروری می‌باشد. هدف از این تحقیق تهیه زیرمدلی 
کاربردی برای تعیین بهترین تاریخ کاشت گیاه کینوا بسته به نوع 


اقلیم و رقم مورد استفاده به لحاظ طول دوره رسیدگی (زودرس؛ 
متوسطرس و دیررس) می‌باشد. 


مواد و روش‌ها 

این آزمایش در دو منطقه از استان یزد (شهر یزد به‌عنوان منطقه 
از سطح دریا ۸۰۰ متر) طبق زمان‌بندی جدول ۱ اجرا شد. آزمایش 
اول در مزرعه شخصی کشاورزی در شهر یزد با طول و عرض 
زراعی ۱۳۹۸-۹۹ اجرا شد و آزمایش دوم در مزرعه شخصی 
کشاورزی در روستای ارنان‌یزد با طول و عرض جغرافیایی ۳۱۳۱۷ و 
شد. نکته قابل توجه در تفاوت اقلیم دو منطقه یعنی اختلاف حداقل 
دماهای این دو ناحیه بین ۸ تا ۱۰ درجه سانتی گراد و اختلاف حداکثر 
دماهای آن‌ها بین ۵ تا ۷ درجه سانتی‌گراد بود (شکل ۱). 
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۱3 2321[ همم عطاصماموه آ۰۹۱۳ 
اسفند بهمن مهر نیمه شهریور. اول شهریور 
29 29 20 6 23 
28 5 12 19 
۲ 
سس ]۷[ ده ده ۳۱01/0101 
0 3 40 
د 35 
0 ِ: 
ِ رز ود 
82 50 ۲ درو و 
ِ 8 20 ه 
٩‏ نو 3 
5 مد | نو ۶ 
-] 4 
20 1 
3 د و ۳ 
10 ۳ 
0 10 


6 0۲ 00۲ ع۹ عسته لد سنا معا( توش معا طاع۳ . و[ 


طاجما( 


(۱۱) صمت)دانورن۱ 


۱/۵ ۸۲1 (7 5 6 ۲ 

مرداد ۳ خرداد اردیبهشت فروردین 10600۵۵ 
۶۴ 216 ۷2270 

26 28 28 29 29 عصتاصمام 


تاریخ کاشت منطقه یزد 
۱ 
عصتاصهام 
تاریخ کاشت منطقه ارنان 


7 سس ۷0( ده ۱0[0[]200ع۳ 


)4 1000۱00۵۵۲۰ 
۱ 
ت 


۶ 0۲ 06۱ ۹6 عنش لگ سنا لتق معا و۲ مقر 


دم 


شکل ۱- خلاصه آمار هواشناسی شهر یزد (۸) و روستای ارنان (1۶) در زمان آزمایش 
سنوی معط ۶ متا عط) 24 (ظ) عمهااز۲۷ مه 0صد ره سا ۷۲۵۸۵ ۵ ماامتاهاه اامماممعام ۵ من .1 مستاهز۲] 


تاريخ کاشت‌های این تحقیق در ۱۰ آزمايش جداگانه برای هر دو 
وه ونور قالت آزمایکی منم بایکهای کال تصاافن باامبه نکر 


انجام شد. عامل‌های آزمایشی شامل ۵ لاين امید بخش اصلاح‌شده 
در مرکز تحقیقات شوری یزد به‌همراه رقم تی‌تی‌کاکا بود. لاین‌ها 


۵۲ نشریه پژوهشهای زراعی ایران. جلد ۱ شماره (۵ ببهار رف 


شامل چهار لاین متوسطرس که در این آزمایش با شماره‌های 
ارطممونن ما۰ م0 و عه00 53لا 
یک لاین دیررس که با شماره ۴ (091001) و رقم زودرس 
تی‌تیکاکا که با شماره ۵ مشخص شده بود. ۳ ماه آبان» آذر و دی در 
یزد و ۴ماه آبان» آذر دی و بهمن در روستای ارنان به‌دلیل سرمای 
زیاد تاریخ کاشتی انجام نشد. از آن‌جایی که تاریخ کاشت اصلی گیاه 
کینوا در منطقه یزده شهریور است (2019 ,.۵1 2 92169 بنابراین 
دو تاریخ کاشت در شهریور (نیمه اول و نیمه دوم شهریور) اجرا شد. 
کاشت به روش خشکه کاری در کرت‌هایی به ابعاد ۲#۲۳متر» با دست. 
به روش خطی و در عمق یک تا دو سانتی‌متری انجام شد و بلافاصله 
بعد از کاشت به‌صورت غرقابی آبیاری شد. بعد از سبزشدن,» تراکم 
بوته به ۶۴ بوته در مترمربع رسید. 

داده‌های هواشناسی از نزدیک‌ترین ایستگاه هواشناسی به محل 
کشت تهیه شد. ایستگاه هواشناسی یزد واقع در فرودگاه شهر یزد با 
طول و عرض جغرافیایی ۳۱۳۵۴ و ۵۳۹۱۶ است و روستای ارنان 
به‌دلیل نداشتن ایستگاه هواشناسی از داده‌های ایستگاه هواشناسی 
روستای گاریز واقع در ۸ کیلومتری این روستا با طول و عرض 
جغرافیایی ۳۱۶۱۸ و ۵۴۹۰۶ استفاده شد. 

نمونهگیری و یادداشت‌برداری‌ها به‌طور مرتب» هر سه روز یک‌بار 
و متناسب با پیشرفت مراحل فنولوژیک هر لاین انجام شسد ( تطعاو8 
9 ,. 67). یادداشت‌برداری مراحل فنولوژیک» مشاهده بوته‌های 
دو ردیف از وسط هر کرت در طول دو متر بود. ملاگ عمل 
یادداشت‌برداری هر مرحله از نمو (سبزشدن, گل‌دهی و تغییر رنگ 


بوته‌ها)» ۵۰ درصد بوته‌های مورد بررسی بود. 

محاسبات آماری با کمک نرم‌افزار 2010 ,12601 انجام شد و 
تمامی نمودارها نیز در محیط 2010 ,آتآ)۲2 ترسیم شدند. 

برای محاسبه درجه-روز-رشد از روابط (۱)» (۲) و (۳) استفاده شد. 


)۱( /(صتص ] + هط ۷۲۳[ 
)۲ 11-11۷۲۳1-60 11 
)۲ 1171-+11-1 ]1-2[ 1 


که در این معادلات 11۷۳ میانگین دمای روزانه» ,م1 حداکثر 
دمای روزانه» مزمم1 حداقل دمای روزانه» 10 دمای پایه» 1211 درجه 
روز رشد روزانه» 1 درجه روز رشد تجمعی» روز مورد نظر برای 
محاسبه روز ذ ام و 1-ذروز قبل از آن هستند ( ,41 6 نصعااه6 
26 

برای به‌دست آوردن پایه اطلاعات درجه-روز-رشد ابتدا درجه- 
روز-رشد محاسبه شده برای لاین‌های کینوا به سه گروه زودرس» 
متوسطرس و دیررس دسته‌بندی شد و میانه داده‌های درجه روز 
رشد هر گروه در جدول قرار داده شد (جدول ۲). پس از محاسبه 
طول‌روزهای ویژه لاین‌های کینواآقدام به تهیه زیرمدل با زیان 
برنامه‌نویسی 751 شد. در تهیه جدول طول‌روز با توجه به روزکوتاه 
بودن کینواء برای تمام لاین‌ها در مدل تا ۱۲/۵ ساعت سرعت نمو 
یک و بعد از ۱۳/۸ ساعت سرعت نمو صفر در نظر گرفته شد. دمای 
پایه در مدل دو درجه سانتی‌گراد در نظر گرفته شد ( ,۵1 6 تلد0ع۳۲ 
2021 


جدول ۲- میانه درجه -روز-رشد لاین‌های کینوا برای نمو کینوا 
عمصنا دمصنیی ]0 20۱۲۵ دوع صهع]۱۷ -2 12۳016 
پارامتر درجه-روز -رشد لاین‌های کینوا 
معا مصتتان ۵۴ رح مععع0 عمط سم 


دیررس متوسطرس زودرس مراحل نمو 
۷ 26[ ادص صین0ع]۱۷ )۵ ات۲۵ ااصمصصم1(6۲۷۵۱0 
"ماهه دوم 7ماهه نخست "ماهه دوم "ماهه نخست "ماهه دوم "ماهه نخست ه‌فاه 
6 0 16[ 6 164 106 6 000 1106 6 164 106 6 000 1106 6 1۳60 106 
۳۱۵0۵06 )۱۵0۵ ۳۱۵06 006 ۳۵0۵5 عط ۳۵ 
925 925 925 925 925 سبزشدن 
0 ۴۱۳0۸۵۲۵۵0۵۵8 
1012٩1 1155.5 1155.5 13338 13338‏ 10012۰1 گل‌دهی 
7 و1وعطاجخ 
رسیدگی فیزیولوژیک 
2357 2357 1061.5 1015.5 1531173 153173 لمعنعهامنووط 
۱/۹۹۹۹ 


با توجه به تفاوت زیاد داده‌های درجه-روز-رشد ۶ ماهه اول سال 


میلادی با شش ماهه دوم برای هر گروه از لاین‌های کینوا ۲ جدول 
درجه-روز-رشد ۶ ماهه اول (از دی تا خرداد) و ۶ ماهه دوم (از تیر تا 


آذر) تشکیل گردید (2021 ,41 نادع). بعد از تهیه زیرمدل نمو 
تهیه شده بارها و در هر بار شش دفعه اجرا شد (سه گروه زودرس» 
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متوسطرس و دیررس و برای هر گروه دو دفمه شش ماهه اول و 
شش ماهه دوم سال) و خروجی آن با اندازه‌گیری‌های شهر یزد 
مقایسه شد. در صورت لزوم تغییراتی در داده‌های جدول درجه-روز- 
رشد مدل داده شد و دوباره مدل اجرا گردید تا بهترین خروجی و 
نزدیک‌ترین داده به مقدار اندازه‌گیری به‌دست آمد که به‌عنوان برآورد 
زیرمدل به زبان 75۲ تبدیل شد (شکل ۲). 

در نهایت ارزیابی زیرمدل با داده‌های روستای ارنان و بر اساس 
پارامترهای برآورد شده در واسنجی انجام شد. 

برای واسنجی و ارزیابی کارایی زیرمدل از جذر میانگین مربعات 
خطا (1۷15۳8)» ضریب تغییرات یا درصد 11517 (46۷ شاخص 
توافق ویلموت (0 کارایی مدل (1۷]۳» مبانگین انحراف از مدل 
(0۷18» ضریب تبیین ( 418 آزمون خط ۱:۱ روز تا سبزشدن» گل‌دهی 
و رسیدگی فیزیولوژیک استفاده شد. 


۵۳ 


جپت مقایسه اختلاف بین مقادیر شبیه‌سازی شده و 
مشاهده شده استفاده شده است که هرچه ۷۲9۳ به صفر نزدیک‌تر 
باشد نشان‌دهنده آن است که شبیه‌سازی بهتر انجام شده است. 
چنانجه مقدار 15718۳ کمتر از ۱۰ درصد باشد شبیه‌سازی عالی» بین 
۰ تا ۲۰ درصد خوب. بین ۲۰ تا ۲۰ درصد متوسط و بیشتر از ۳۰ 
درصد ضعیف ارزیابی خواهد شد (2019 ,.0 61 255171). هرچه 
مقدار 4 به‌دست آمده توسط مدل به یک نزدیک‌تر باشد بیانگر آن 
است که مدل در شبیه‌سازی موفق‌تر بوده و توانسته است میزان 
پیش‌بینی صفت را با اختلاف کمتری نسبت به مقلدیر مشاهده شده 
انجام دهد (1982 ,۱۷111011). بهترین حالت 1۷8 صفر است و 
هرچه از صفر کمتر باشد یعنی مدل کمتر از وافعیت و هرچه از صفر 
بزرگ‌تر باشد یعنی مدل بیشتر از واقعیت تخمین می‌زند. شاخص 
کارایی مدل 2۷۲۳ بین صفر و یک خوب است و بهترین حالت عدد 
یک است و بدترین حالت؛ زمانی که منفی باشد. 
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شکل ۲- دستورات زیرمدل نمو کینوا به زبان برنامه‌نویسی 51 بر اساس درجه-روز -رشد 
(61۱) ۲۵ مرمع همز 60۵ رن موه ءمهممه۱ ۲٩۲‏ ضز م0صمصصصصی عم امعم مام م0 ممصتام) -2 موز ۲ 


و( نشریه پژوهشهای زراعی ایران. جلد ۱ شماره [3 ببهار ۱:۰۲ 


ضرایب رگرسیون ساده خطی (۸ و 0) بین مقادیر پیش‌بینی شده و 
مقادیر واقعی» ضرایب 2 و 9 به‌ترتیب نشان‌دهنده مقدار انحراف خط 
رگرسیون از مبداً مختصات و مقدار شیب خط رگرسیون از خط ۱:۱ 
نشان‌دهنده ایده‌آل بودن مدل است. خط ۱:۱ دارای عرض از مبداً 
صفر (220) و شیب ۴۵ درجه (0-1) است (2019 ,.۵1 6 تتلهع2ل). 


زیرمدل تهیه شده به زبان 75 برای پیش‌بینی سه مرحله اصلی 
نمو لاین‌های کینوا (سبزشدن, گل‌دهی, رسیدگی فیزیولوژیک) در سه 


گروه زودرس» متوسطرس و دیررس آماده شد. ابتدا پارامترهای 
زیرمدل با استفاده از میانه داده‌های درجه-روز-رشد واقعی مراحل 
فنولوژیک به‌دست آمده در آزمایش‌های مزرعه‌ای شهر یزد استخراخ 
شدند (جدول ۲ سپس عمل واسنجی یعنی تنظیم پارامترهای جدول 
درجه-روز-رشد زیرمدل در محدوده معقول آن‌ها انجام شد به 
طوری که نتایج شبیه‌سازی شده, قابل مقایسه با مقادیر اندازه‌گیری 
شده باشد. در نهایت داده‌های جدول ۲ به شرح جدول ۳ تغییر يافت 
باتوی که یبش فراع اصلن تقو با عناق مقباز شناکض هاش 
سنجش آماری صورت گیرد. 


جدول ۳- مقادیر درجه-روز-رشد لاین‌های کینوا بعد از واسنجی زیرمدل نمو کینوا 
1 همین 0 ممتاهدانای ماه ممصنا مین 0 ممتاه (6۵) رون مرمع( مز06۳0۵۳۷ -3 12۳016 
پارامتر درجه-روز -رشد لاین‌های کینوا 
ععصنا ومصشای 0۴ رح م0626 عصنطا م6 


دیررس متوسط رس زودرس مراحل نمو 
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ننایج اجرای زیرمدل تهیه شده به زبان ۳5۲ و مقدار شاخص‌های 
آماری برای تعداد روز از کاشت تا گل‌دهی در شکل ۳ ارائه شده است. 
برای تعداد روز از کاشت تا گل‌دهی بر اساس این ننایج» مقدار جذر 
میانگین مربعات خطا (۷]5۳) در دامنه ۵/۵ تا ۶/۵ روزه ضریب 
تغییرات یا درصد ۴۷۲۹۴ (0۷) در دامنه ۹/۵ تا ۱۳/۵ درصد. 
شاخص توافق ویلموت (1) در دامنه ۰/۸ ۱ ۰/۹ کارایی مدل (۷۳) 
در دامنه ۰/۴۳ ۱ ۰/۷۶ میانگین انحراف از مدل (۷۲8) در دامنه 
۲۹ و ضریب تبیین (8) در دامنه ۰/۸۳ تا ۰/۹۲ قرار 
گرفت. بنابراین مرحله گل‌دهی بسیار خوب تخمین زده شده است. 
نتایج آزمون ) مربوط به مقایسه خط رگرسیون برازش داده شده با خط 
نیز نشان داد که شیب و عرض از مبدا خط رگرسیون با خط ۱:۱ 
اختلاف معنی‌داری نداشت (شکل ۳). 

در بررسی تعداد روز از کاشت تا گل‌دهی» بزرگ‌ترین ضریب تبیین 
معادله رگرسیون در لاین‌های دیررس (۹۳/--82) و کوچک‌ترین 
ضریب تبیین در رقم زودرس (۰/۸۶-*1۳) بهدست آمد. این امر نشان 


از میزان پراکندگی نقاط در شبیه‌سازی این مراحل نموی دارد و بدین 
معنا است که با معادله رگرسیون ساده خطی در رقم زودرس ۸۶ درصد 
و در لاین‌های دیررس ٩۳‏ درصد از تغییرات در مقادیر شبیه‌سازی‌شده 
تعداد روز از کاشت تا گل‌دهی, به‌وسیله تغییرات در مقادیر مشاهده‌شده 
تعداد روز از کاشت تا گل‌دهی کینوا قابل توضیح است (شکل ۲). 

براساس شاخص میانگین انحراف از مدل (۷189)» تطابق 
داده‌های شبیه‌سازی‌شده روز تا گل‌دهی با مقادیر مشاهده‌شده در لاین 
دیررس کینوا بهتر از سایر لاین‌ها اجرا شد (شکل ۲). 

در بررسی شاخص کارایی مدل (۷1۳ در تمامی موارد مثبت شد 
بنابراین مدل در شبیه‌سازی روز تا گل‌دهی قابل‌قبول است. بهترین 
حالت شاخص کارایی مدل (1۷177-۰/۷۶) در لاین‌های متوسطرس 
محاسبه شد زیرا در این شاخص هرچه عدد به یک نزدیک‌تر باشد 
خطای شبیه‌سازی» کمتر است (شکل ۳). 

بالاترین مقدار شاخص توافق ویلموت (*/2۰) در رقم زودرس 
تی‌تی‌کاکا و کمترین مقدار شاخص توافق ویلموت (1<۰/۸) در 


اقبالی شاه‌آباد و همکاران. شبیه‌سازی نمو فنولوژیک کینوا بر اساس داده‌های مزرعه‌ای ۰ ۵۵ 


لاین‌های متوسطرس محاسبه شد که نشان از شبیه‌سازی خوب تحقیقی دیگر نتایج واسنجی تعداد روز تا گل‌دهی» توسط پراگر و 
داده‌های روز تا گل‌دهی» در تمامی لاین‌های کینوا است (شکل ۳ مارا تام یرای وف شم سین رنه 

نتایج واسنجی تعداد روز تا گل‌دهی با مدل‌های تجربی در دماو کینوا در آلمان نشان دادند که مقدار جذر میانگین مربعات خطا 
طول‌روزهای متفاوت برای کینا نشان داد که مقدار جذر میانگین ‏ ۰ (020058) ۳/۱ روز و شاخص توافق ویلسوت (0) ۰/۹۸ برای هر دو 
مربعات خطا (1۷]5۳) ۷/۲ تا ۷/۸ روز و ضریب تبیین (*) ۰/۷۹ تا رقم به‌دست آوردند که به نتایج این تحقیق نزدیک است ( 6 ۳۳۵۵6۲ 
۵ به‌دست آمد. در بررسی شیب و عرض از مبدا خط رگرسیون با 09 ,.01). 


خط ۱:۱ اختلاف معنی‌داری نداشت (19992 ,.]» 61 6۲06۲0 در 
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شکل ۳- مقادیر مشاهده‌شده و شبیه‌سازی‌شده روز تا گل‌دهی منطقه یزد در رقم زودرس (۸)؛ لاین‌های متوسطرس ()» لاین دیررس (6) و 
لاین‌های کینوا (10) به همراه جذر میانگین مربعات خطا (418۱1517 ضریب تغییرات یا درصد 038۱1518 (60۷), شاخص تواف ویلموت ()» کارایی 
مدل (۷110)؛ میانگین انحراف از مدل (0119)» ضریب تبیین (*18) خط 1:1 (صاف) و خط ر گرسیون (بریده) 
و(() کعصنا )تحص تلع ررض تایه امه ما ممنعع۲ ۷۵2۵ ۵۲ ممته10۳۲۷ و) مرج ۵ ما۲ 60اها نی مره 0ع۲مووا) -3 مساهز۲] 
و(62۷) ۱۱۲۸۱۷۲۹۲۶ ۷۵۱۵۸10 ۵۲ )صمت‌تالهی و(تاها1۵۷) مه هنود همه )م۳۵ (ظ) معا مصتمن 20 (0) مصنا راتسامصه ما12 


1 10۵ ,22) «صصنصهع0 ۵۴ اممت‌تالهی و(ظ۷۲۳) ممتاه06۷1 عم همم ور مت له مه ررم) 06ص غصمصصمعمه )ممصاز ۱۷۷ 
(اجه) مصتا رمتقعع۵۱ع۳ 0صد (طاممصصم) 


به‌منظور ارزیابی و اطمینان از کارایی زیرمدل تهیه شده داده‌های جذر میانگین مربعات خطا ( 5۷]5۳) بین ۴/۴ تا ۶/۴ روز شاخص 
روز تا گل‌دهی شبیه‌سازی شده با داده‌های مشاهده شده آن‌ها در توافق ویلموت (1) بین ۰/۹۹ تا ۰۱ کارایی مدل (۷]۳) بین ۰/۹۶ تا 


ضریب تغییرات یا درصد ۷۲5۳ (۷)) بین ۷ تا ۱۲ درصد. تبیین (*3) بین ۰/۹۲ تا ۰/۹۸ بود. این مقادیر نشان از برآورد خوب 


2۶ نشریه پژوهشهای زراعی ایران» جلد ۱ شماره (۵ ببهار رف 


روز تا گل‌دهی کینوا با زیرمدل نوشته شده به زبان ۳۹۲ داشت و 
مقادیر روز تا گل‌دهی شبیه‌سازی شده اعتبار لازم را کسب کردند 
(شکل ۴ 

مقایسه شیب خط رگرسیون بین مقادیر مشاهده‌شده و 
پیش‌بینی‌شده با خط ۱:۱ در مورد روز تا گل‌دهی کینوا نشان داد که 
اختلاف معنی‌داری بین شیب خط رگرسیون و شیب خط ۱:۱ وجود 
نداشت؛ بنابراین زیرمدل واسنجی‌شده از اعتبار کافی برای شبیه‌سازی 
روز تا گل‌دهی کینوا برخوردار بود. 

نتایج ارزیابی تعداد روز تا گل‌دهی با مدل‌های تجربی در دماو 
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طول‌روزهای متفاوت برای کینوا نشان داد که مقدار جذر میانگین 
مربعات خطا (80۷]5۴) ۴/۵ تا ۴/۹ روز و ضریب تبیین (*) ۰/۸۵ تا 
۷ به‌دست آمد. در بررسی شیب و عرض از مبدا خط رگرسیون با 
خط ۱:۱ اختلاف معنی‌داری نداشت (19992 .۵ 1 71070ظ). 
محققین دیگر نتایج ارزیابی تعداد روز تا گل‌دهی بامدل 05۲ 
۵0 برای دو رقم جسی و زنو کینوا در آلمان نشان دادند 
که مقدار جذر میانگین مربعات خطا (11۱۸515) ۱/۵ روز و شاخص 
توفق ویلموت (0) ۰/۹۹۷ برای هر دو رقم بهعدست آوردند که به نیج 
این تحقیق نزدیک است (2019 ,.4 61 ۳۳2826۲ 
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شکل - مقادیر مساهده‌شده و شبیه‌سازی‌سده روز تا گل‌دهی منطقه ارنان در رقم زودرس (۸۵)ء لاین‌های متوسطرس (1), لاين دیررس () و 
لاین‌های کینوا (10) به همراه جذر میانگین مربعات خطا (418۱1517 ضریب تغییرات یا درصد 018۱1518 (60۷), شاخص تواف ویلموت ()» کارایی 
مدل (۷11» میانگین انحراف از مدل (40۷119 ضریب تبیین (*418 خط ۱:۱ (صاف) و خط رگرسیون (بریده) 
کعصن تاحفص تلع رره) نها تایه «ابهه ص ممنعع۲ محصم ۵ مصته۳7 10 ما رل ۵ ماه 60) هی 0صره 0و0 -4 مساهز۲] 
,(0۷) ۱۱۱۸۱۷۲۹۲ م)ه۱ه۲ ۵۴ )مهن و(1۱۷]۹۳) مه مهن وی صممجه عممنا از (ظ) ممصتا هم 0صو (0) مصنا راتساطد ماها ررظ) 


1 :1 ,12) «صصنصهع0 ۵1 اممت‌تالهی و(ظ۷۲۳) «متاه06۷1 عم همم ور تال مه ررم) 06ص غصمصصمعمه )ممطز ۱۷۷ 
(اجه) مصتا متقعع ۳6۵ 0صد (طاممصصم) 


مقایسه نتایج زیرمدل با داده‌های اندازه‌گیری شده و مقدار 
قلخضی‌های آماری برای تاذ روز تا رسیدگن فیزیولیژیک دز شکن ۵ 
یسک بان تاه روا زگ ریب نک بر آنانن آنتن 
نتایج مقدار جذر میانگین مربعات خطا (18۷۹۳) در دامنه ۱۰/۶ تا 
۹ روزء ضریب تغییرات يا درصد ۱۴۷۲۹۴ (0۷) در دامنه ۱۰/۵ 


تا ۱۴ درصد» شاخص توافق ویلموت (0) در دامنه ۰/۶۶ ۰۰/۶۸۱ 
میانگین انحراف از مدل (۷]8) در دامنه ۰/۰۹ تا ۰/۱۱و ضریب تبیین 
() در دامنه ۰/۸۸ تا ۰/۹۱ قرار گرفت؛ بنابراین مرحله رسیدگی 
فیزیولوژیک بسیار خوب تخمین زده شده است. نتایج آزمون ] مربوط 
به مقایسه خط رگرسیون برازش داده شده با خط ۱:۱ نیز نشان داد که 


اقبالی شاه‌آباد و همکاران. شبیه‌سازی نمو فنولوژیک کینوا بر اساس داده‌های مزرعه‌ای ۰ ۵۷ 


شیب و عرض از مبدا خط رگرسیون برازشی با خط ۱:۱ اختلاف 
معنی‌داری نداشت (شکل ۸۵). 

در بررسی تعداد روز از کاشت تا رسیدگی فیزیولوژیک» بزرگ‌ترین 
ضریب تبیین معادله رگرسیون در لاین‌های دیررس (۳-۰/۹۱) و 
کوچک‌ترین ضریب تبیین در لاین‌های متوسطرس (۰/۸۸-*18) 
به‌دست آمد. این امر نشان از میزان پراکندگی نقاط در شبیه‌سازی این 
مراحل نموی دارد و بدین معنا است که با معادله رگرسیون ساده خطی 
در لاین‌های متوسطرس ۸۸ درصد و در لاین‌های دیررس ٩۱‏ درصد 
از تغییرات در مقادیر شبیه‌سازی شده تعداد روز از کاشت تارسیدگی 
فیزیولوژیک به‌وسیله تغییرات در مقادیر مشاهده‌شده تعداد روز از 
کاشت تا رسیدگی فیزیولوژیک کیناقبل توضیح است (شکل 4۵ 

بر اساس شاخص میانگین انحراف از مدل (۷۲5) تطابق داده‌های 
شبیه‌سازی‌شده روز تا رسیدگی فیزیولوژیک با مقادیر مشاهده شده در 
لاين دیررس کینوا بهتر از سایر لاین‌ها اجرا شد (شکل ۵). 

بالاترین مقدار شاخص توافق ویلموت (02۰/۶۸) در رقم زودرس 
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تی‌تی کاکا و کمترین مقدار شاخص توافق ویلموت (12۰/۶۶) در 
لاین‌های متوسطرس محاسبه شد که نشان از شبیه‌سازی خوب 
داده‌های روز تا رسیدگی فیزیولوژیک در تمامی لاین‌های کینوا است 
(شکل ۵). 

نتایج واسنجی تعداد روز تا رسیدگی فیزیولوژیک پراگر و همکاران 
با مدل 030۳00-](05 برای دو رقم جسی و زنو کینوا در آلمان 
نشان دادند که مقدار جذر میانگین مربعات خطا (۴(۷]۹۳) ۷ روز و 
شاخص توافق ویلموت (8) ۰/۸۹۵ برای هر دو رقم به‌دست آوردند که 
به نتایج این تحقیق نزدیک است. پیش‌بینی نمو فنولوژیک يا مراحل 
نموی گیاهان اهمیت زیادی دارد چون تولید و توزیع ماده خشک در 
مدل‌های شبیه‌سازی گیاهان زراعی تا حدود زیادی تحت تاثیر زمان 
وقوع مراحل فنولوژی می باشد. به همین ترتیب, ممکن است 
پیش بینی تغییرات سطح برگ و ساير فرآیندها به پیش‌بینی نمو 
فنولوژیک وابسته باشند (2006 ,۵1 61 نصهاا60). 
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شکل ۵- مقادیر مساهده‌شده و شبیه‌سازی‌شده روز تا رسیدگی فیزیولوژیک منطقه یزد در رقم زودرس (۸/؛ لاین‌های متوسطرس (؛ لاين 
دیررس ()) و لاین‌های کینوا (1) به همراه جذر میانگین مربعات خطا (181517)» ضریب تغییرات یا درصد 18۷1517 (۷))» شاخص توافق ویلموت 
(0)» کارایی مدل (110)؛ میانگین انحراف از مدل (41118 ضریب تبیین (*408 خط ۱:۱ (صاف) و خط رگرسیون (بریده) 
نع رره) وتان امه ص ممنمع۲ ۷۲۵۸۵ ۵۶ تسام ام‌تومامتعودام ما درل ۵۶ ممتاه۲ 0م)هاصته 0ج ۲۵0 0و0 -5 عاع۲1 
0600 و(۱۱۷۲۹۲) مه موز عم متا از () ممصنا میتی 0صج ر)) مصن راتستااحصه عاها ررظ) عمصنتا پاتست/م 


0 )هم و(0۷۲) 06۲۷۱2600 ۵061 همه وا مت که ۵06 ور) 1۱06 )صمصرعع هه )مصصل ۲۷۲۷ ,(0۱۲۷) ,۱۸۱۷۲۹۲ 
(اسه) عصنا جمنویی ۲6۵ 0صد (طاممصع) 1:1 عصنا ر۳) 


2۸ نشریه پژوهشهای زراعی ایران» جلد ۱ شماره (۵ ببهار ۱:۰۲ 


به‌منظور ارزیابی و اطمینان از کارایی زیرمدل تهیه شده, داده‌های 
روز تا رسیدگی فیزیولوژیک شبیه‌سازی شده با داده‌های مشاهده‌شده 
آن‌ها در آزمایش میدانی منطقه ۲ (روستای ارنان) مورد ارزیابی قرار 
گرفت. 

ضریب تغییرات یا درصد 0۴۷۲5۳ (۷)) بین ۶/۸ تا ۸/۶ درصد. 
جذر میانگین مربعات خطا (1151) بین ۶/۲ 1 ۸/۷ روز شاخص 
توافق ویلموت (8) بین ۰/۷۵ تا ۰/۹۲» میانگین انحراف از مدل (2158) 
بین ۰/۰۵ تا ۰/۰۸ و ضریب تبیین (**) بین ۰/۹۲ تا ۰/۹۸ بود. این 
مقادیر نشان از برآورد خوب روز تا رسیدگی فیزیولوژیک با زیرمدل 
نوشته‌شده به زبان ۲۹۲ داشت و مقادیر روز تا رسیدگی فیزیولوژیک 
شبیه‌سازی شده اعتبار لازم را کسب کردند (شکل ۶). 
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مقایسه شیب خط رگرسیون بین مقادیر مشاهده‌شده و 
پیش‌بینی‌شده با خط ۱:۱ در مورد روز تا رسیدگی فیزیولوژیک کینوا 
نشان داد که اختلاف معنی‌داری بین شیب خط رگرسیون و شیب خط 
۱ وجود نداشت؛ بنابراین زیرمدل واسنجی‌شده از اعتبار کافی برای 
شییه‌سازی روز تا رسیدگی فیزیولوژیک کینوا برخوردار بود. 

نتایج ارزیابی تعداد روز تا رسیدگی فیزیولوژیک سایر محققین با 
مدل 053-070۳0180 برای دو رقم جسی و زنو کینوا در آلمان 
نشان دادند که مقدار جذر میانگین مربعات خطا (1577) ۷ روز و 
شاخص توافق ویلموت (8) ۰/۸۹۵ برای هر دو رقم به‌دست آوردند که 
به نتایج این تحقیق نزدیک است (2019 ,.اه 6 ۲۲2826۲ . 
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شکل "- مقادیر مشاهده‌شده و شبیه‌سازی‌شده روزها تا رسیدگی فیزیولوژیک منطقه ارنان در رقم زودرس ۳۹ لاین‌های متوسط رس (ظ)» لاین 
دیررس () و لاین‌های کینوا () به همراه جذر میانگین مربعات خطا (18151) ضریب تغییرات یا درصد 18۷151 (0۷), شاخص توافق ویلموت 
(0)» کارایی مدل (110)» میانگین انحراف از مدل (41118 ضریب تبیین (4087 خط ۱:۱ (صاف) و خط رگرسیون (بریده) 
منم رزه) تایه واه صز ممتعع فص ۵ اتمه آی‌تعماملوودام آتاصس مرمل اه ممتاه۲ 60) هی 0صره ۲0 موها0 6 عستاه۲1] 
۲۵۰۱۵0۱ 0۴ انامه و(۱۷]۹۳) مه موی همه اممنا ۱۷۱۵۵ (ظ1) مم‌صنا مصتان 20 (0) 6ص مها ررظ) معصتا ۸2۲0016 


مصنا ,02 صصنتصماع0 ۵ )نی و(ظ۷۱) 0612610 ۵061 همه ور متا ۳۵061 ورم) 106 امه عمصصا ۲۷۲۷ (0۱۷) 
(اجه) 6ص1 ممنوعع 6۵۱ 290 (رطاممصصرع) 1 :1 


اقبالی شاه‌آباد و همکاران. شبیه‌سازی نمو فنولوژیک کینوا بر اساس داده‌های مزرعه‌ای ۰ ۵٩‏ 


از آن‌جایی که درجه-روز-رشد. پارامتر متداول و کاربردی در 
پیش‌بینی مراحل مختلف فنولوژی گیاهان است (2000 ,255171]) به 
کمک آن می‌توان تاریخ کاشت مناسب را مشخص نمود و در ساخت 
مدل‌های مکانستیکی ساده قابل استفاده است. در آزمون‌های انجام 
شده توسط محققین مشخص شد که قدرت مدل‌های ساده بیشتر از 
مدل‌های پیچیده است (2021 ,اه ۶ «92ط2). مدل‌های 
مکانستیکی ساده علاوه بر این که محدودیت مدل‌های تجربی به 
لحاظ کاربرد در اقلیم‌های مختلف نداشته و به هزینه‌های زیاد برای 
تعیین پارامترهای مناطق جدید هم نیاز ندارند. در ضمن تهیه 
واسنجی و اعتبارسنجی زیر مدل مکانیستیکی نمو کینوا کمکی در 
جهت ساخت مدل کامل شبیه‌سازی رشد کینوا در آینده می‌باشد. 


نتیجه گیری 


قا مایخ مطالیت نام واه که ملال نو تیه تم یزاین که بت 
حسب درجه-روز-رشد به خوبی مراحل نمو (سبز شدن» گل‌دهی و 
رسیدگی) انواع این گیاه به لحاظ رسیدگی (زودرس, متوسطرس و 
دیررس) پیش‌بینی می‌کند و می‌توان به کمک آن تاریخ مناسب 
کاشت هر متاطی معا شین کرد ازنن زیرمتال کتالبیره شوه 
هم‌اکنون می‌تواند برای ارزیابی تاثیرات مختلف دمایی و فتوپریودی 
برای تصمیم‌گیری در طیف وسیعی از محیطهای رشد در نظام‌های 
کشت کینوا در شرایط اقلیمی جاری و آینده مورد استفاده قرار گیرد؛ 
بنابراین از این زیرمدل در کارهای آموزشی- تحقیقاتی و کاربردی در 


مزرعه می‌توان استفاده کرد. 


ز 3 ۱۰ 


۵ 6۱۵۵0۴۵001۱۱ ۱۵۷ 0۴ ۳2۵ 8616116 0۴ ۷۵/۷۵/۱۵۷۱ 2018(۰) ,۱۷۲ باتعطعظ ع بخ طلمجصلته رو رلفدطاظ 
ععصمنهک عععتهمم صمتصم! لمممناد]( 15۳ .عنم آمعومامام موه وه عمط دمبایی (۲۷۷۵/۱۵۰ موی نکوم0۵90) 
۵0۰ طوتام م۴ طاز۳ صفتویعظ صل) بصق رزح‌تهک :2019 ۹6۵۰ 2-5 رووعتهومن) 

مه 1۵0۷ ۵۶ ووفممموع۲ تمصع مه 0۵۵100)مطه عصذ1[ع۷]۵0 .(19992) .[ .۸ بالق ,۷۷۰ بک رلک .لا .1 ,۳۲6۲0 
19-4۰ ,63 ,6۶۵۲6۲ ۵۵ ۳۲۶۵/۵0 ۷۷۱۱۱۸۰(۰ 0۱۱۵۵ 0۷۱۵۵۴۵00) دمصننان ظ1 

همصنیان صا فعفقطام ممممماهبه 910۷۵همو-م۵مم)مطظ ,(19990) .۲ یه ملق عک .۷۷ نک ,لک .لا بط ,۲6۲0ظ 
231-۰ ,60 ۵۹۵۵۲۵1 ۲۵۵۶ ۳۸۵۱۵4 ۱۷۷۱۱۱۵۰(۰ ۱۱۵۵و 0۵00 0۵)) 

هممننان گه ومعهاد لممنعماممعطم ۵۶ ععممموعع 1 .(2021) ,۷۲ ,تطماهه ع ,۷ ,۱۷۲ رتتلعع۱2 ,۷ مصقطعل یط بئله‌ططعظ 
61-۰ ,9 7 ,660۲6 ۲۵۵5 ۲۱۵۱ ۵۴ ۵۱۱۵ ۲۳۵/۱۱۵۱ .عمصطنعع۲ ماما 0ص وعناجتهم‌جصع) م0 فعصناً عصوتططم۳0م 
2 101.01271026 //:0009 

27 ۸۵۵62560 موم 00و۵۳ 10۲۵و که ۲۵۵0 ۱۱۵۱۵۵۵۲5۰ 171 ۵۵۲۵0 ۵۳00 ۳۲۵۵۵ 2014(۰) ۲۸۵۵۰ 
20015۰ 

)۲۵ ,(2017) .1 ,0921 ی .۲ رلق1ت وبا معاصقصصه وا ظ .۷۷ مق210171 ۷ و۷ مقصهت و.ل م12 ویک مقطع1101 ۷۷ .9 ,۲1506۲ 
٩۲655‏ ۷۵۵۲ ۱0 0ماانصهاتاه ۵۵05و حصم انا۵۲0٩‏ (۷۷۱۱۱۵۰ ممیی 0۱۵۱۵0۵۵) دمصتیان ۵۶ ممتالدممجمع )صهدمتاصه هه 
.558-4 ,107 ۳۳0۵ ك ۳905 ۱/۵ هیهت م11 200 «اتصنلمه 
005://001.012/10.1016/[.10:06700.2017,/(5 

امعصمماه 0۵7 مه وعاهل م90۱۳ 01 ]۱۱۵۲۵۵۵0-۳۲]6۵ صا مصلنام 2014(۰) .ظ .9 رجع‌وطا0ع12 عک ریک مطقا0۱10-۵۱ظ) ربظ بطعت۲۱ 
1 8 ۸105://001.018/10 .1-1 ,23 ,0206 ۳۵۵ ۵0 موه 0 ۵۵۱۱۵ ۷۱۵۱0۰ 20 

40ج دوم عصتاوم۸0 ,(2013) .م1 .ظ بااعاه۳۱ ع ,ین بخ بصماوا نگ هط با رجماعصدنط ر.ظ یک رفسمه 
2 566۳ ۲۵۵0 6۵۵۱ .عصتلاه ۱۵0 ومع ۵ ۲۵1۵ عمط :ات56۵ 1۵00 عتباعمه ما ومتقصصتاه مسابظ ما فصماولاه مصتمم0ه 
89 ,01009://001.012/10 .24-28 

۰ ۱۷۲۸۵320 تطمع طوم‌صدنا مطه[ . وععوعع۳۵م ۶سمع مومت 50۱67۵1 ۸0۵6/۱۱8۵ 2000(۰) ۱۰ ,۷۲ ,۸951۲1 


مج ۳۵۲۵۱۱2۵۲ مموم۱ ۵۶ ممتام2تصت0ام0 ,(2019) .ظ ۸۲ رتتتض ع و رتتاهمط۳2(12 ریک ب بتاعط۵۵0 رب ,۷۲ متتلوعهل . 


بآ ,۸8۵۲0۵60108 0۴ ۱۱۵ ول .2 لوعن۱ مالوممجمن امن ها ممه۷ ۵0 ( اععه س1۳11) ۷۷۳۲۵۵۲ م1 صمتاهع۳۳1] 
72 2 .515-530 


عمجم ۵1 اتاهتام) چم ممقاظ صه مصطا ]1 عصتاصه۱ ۵0۲ ۲۲۲۵۵۲ 2018(۰) ,۷۲ رطمهه2نصوه۱۱ ک وب بتل۷۵12 ,لا رطفت‌ممطونام۳ .. 


679-۰ ,20 ۱0۳۵۳۵۷۵۵4 کم۳0 . صقططریم . رتلفلتدظ .ها عمتاعهبتاانت .10۲ . وعتامنج ۷ .»ما9 صم .. )صصصمترظ 
 ( 9‏ ۵( 
۶ وعووم۵ع0۲۵ 06۷۵۱0۵8۵۵ هه طامع عصتاهلنه ,(2019) .0 .و رعومتمما۲ عک ویو رمصت ول یک ,عا00ظ بیظ ,۳۲۸261 
9 ,۸۵۲0۲0۲۲۲ .۱۱۵061 کی-9 عمط ۵۶ ممتاهبطدبه مه مماهام20 ۱۷۷۲۱۱۱۵۰(۰ ممویی سخمم۵۵)) دمصتبام) 
2 ۵ .832 

0 ۶ ۷۱۵۱0 مه طامع ۵۶ مموتنهم‌جصمن 2003(۰) ۱۷۲ رهم۷۵02 عک بلاط متتاقا نک متصهلمط 0۵ رنه ما0۵ 
1-7۰ ,(3)1 ,6509۲6۵5 ۷۵۵۲۵ ۵0۳۱0 ۵62 ۸۵۱6۵/۱۵۵ 6۵۲9۵۱ .عصمتاتلجمی 160 صنه؟ صا دعتامته۷ 

40 ,۷۱۵10 ,۲۵۵۵۸۵۱۵2۷ ۳۲۵0۱6۲ ۵ م۷۲۵0 9۷10و عصلونا 2021(۰) .۲ متلفطام 2‏ .و بلطهاا0و ر.ظ م1010 ,.9 مصعه‌اطف] 
15-۰ رد3 ۸8۲۵۵۵۵۵8 ۵۴ ۵و جمتانهه) صوع1 صا رما امه معتععهبظ) قامصهت ۵ رهظ بماه ۱۷۷ 
7 ( 


۶ توص ۵۴ 5۵۵0۲ ۲۵۹۵۵۲6۵۵ ۷۵۲۵۲ ]5۵ ۵۲ 66۲۵۵ ی اک ,۵۵6 2017(۰) ۲۰ ,تصح‌طاعطعن1 یک ,۷۲ رتطعل52 
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15۰. 


۶.۰ نشریه پژوهشهای زراعی ایران» جلد ۱ شماره [۵ ببهار رف 


۰ 0۵ تاه مخ 

۶ ۷۱۵۱0 20 وععهاه لهعلعمامصمطام چم ع2 عصتاصهام ۶ه 11601 1 .(2019) .۲ ,رصهطاوده(۱ عک .۷ متصقااهو ,۷ رتطعله 
923-۰ .و2 ,96۵06۵8 ۲0۵ .1 کفعوعه .۲۵۱۵۵ .عصمنانلجهع... فصنلهه ...ها 08ععو .. همصتتان 
1 005://001.0182/10.22077/6868.2019 

۱0۵62م ۱۷۲۱00۵112 2006(۰) بخ متکله۵۵-1>212 ی .۲ ب24ز24-۲6صصهطم] ۷ ر.ل ,۱۷ مط0و)ت06م بیظ ملصه)501 
.14-3۰ ,99 . 1656۲6 . ۲۵88 .. ۳618 .۰ .380860800806 ...م24 صمتام‌المم ...162۴ .. :اصمصصمم781ع0 ...۰ 0هه 
5 ۰ (005://001.012/10.1016/1 

فص دامج؟ ممتاتطت ۲ .(2010) به .ظ رمعطات۱۷]2 ۶ وبا رماصعمناط ر.ظ ب106ل] و.ل ۷۵۲۵۵۵ م۷ بقله۱۷]12 ببظ ,۱۷۵۵2-01۷62 
۴ 0۱۵ 17:6 ۲6۷1۵۷۷۰ 2 :صلدتع صوعصض صع1عصح حح ,(۷۷۱۱۱۵۰ م۵ایای 0۱۵۵۵۵0)) همصتان ۵۶ احتاصماوم لمع0)10] 
8 01/10 .2541-2547 ,90 ,661166 ۸۵۲۱6۵۲۵ 

0 ۱۱۶ 0۴ ۱0/۱29۱ظ .عمصمصصزموم موه گم حمتامتلهبه فطع مم. فاصمعصصمصمع معصمو ,(1982) .۲ .۲ بتامصااز/۷ 
4 ۷ .13 1309-13 ,63 ,»0 ۸۲۵۱20۲۵01081661 

۶ وموطوممووم آهمزو۳۵۹۱0۱0 مه م0 2016(,۰) .ظ .9 بطموها0ع12 عک وی بلهطا10 ۷ بلهژظ بو .9 متقاطاظ بط م۷218 
45-۰ و202 رععععص م0 که 00و 0 ماو .8وعا9 . موجه[ هه تطونهعنا. هه همصتبام0 
3,7 ( .10 /1005://001.0182 
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